Statistical non-parametric mapping in sensor space

Establishing the significance of observed effects is a preliminary requirement for any meaningful interpretation of clinical and experimental Electroencephalography or Magnetoencephalography (MEG) data. We propose a method to evaluate significance on the level of sensors whilst retaining full temporal or spectral resolution. Input data are multiple realizations of sensor data. In this context, multiple realizations may be the individual epochs obtained in an evoked-response experiment, or group study data, possibly averaged within subject and event type, or spontaneous events such as spikes of different types. In this contribution, we apply Statistical non-Parametric Mapping (SnPM) to MEG sensor data. SnPM is a non-parametric permutation or randomization test that is assumption-free regarding distributional properties of the underlying data. The method, referred to as Maps SnPM, is demonstrated using MEG data from an auditory mismatch negativity paradigm with one frequent and two rare stimuli and validated by comparison with Topographic Analysis of Variance (TANOVA). The result is a time- or frequency-resolved breakdown of sensors that show consistent activity within and/or differ significantly between event or spike types. TANOVA and Maps SnPM were applied to the individual epochs obtained in an evoked-response experiment. The TANOVA analysis established data plausibility and identified latencies-of-interest for further analysis. Maps SnPM, in addition to the above, identified sensors of significantly different activity between stimulus types.

MEG and EEG dipole clusters from extended cortical sources

Data from magnetoencephalography (MEG) and electroencephalography (EEG) suffer from a rather limited signal-to-noise-ratio (SNR) due to cortical background activities and other artifacts. In order to study the effect of the SNR on the size and distribution of dipole clusters reconstructed from interictal epileptic spikes, we performed simulations using realistically shaped volume conductor models and extended cortical sources with different sensor configurations. Head models and cortical surfaces were derived from an averaged magnetic resonance image dataset (Montreal Neurological Institute). Extended sources were simulated by spherical patches with Gaussian current distributions on the folded cortical surface. Different patch sizes were used to investigate cancellation effects from opposing walls of sulcal foldings and to estimate corresponding changes in MEG and EEG sensitivity distributions. Finally, white noise was added to the simulated fields and equivalent current dipole reconstructions were performed to determine size and shape of the resulting dipole clusters. Neuronal currents are oriented perpendicular to the local cortical surface and show cancellation effects of source components on opposing sulcal walls. Since these mostly tangential aspects from large cortical patches cancel out, large extended sources exhibit more radial components in the head geometry. This effect has a larger impact on MEG data as compared to EEG, because in a spherical head model radial currents do not yield any magnetic field. Confidence volumes of single reconstructed dipoles from simulated data at different SNRs show a good correlation with the extension of clusters from repeated dipole reconstructions. Size and shape of dipole clusters reconstructed from extended cortical sources do not only depend on spike and timepoint selection, but also strongly on the SNR of the measured interictal MEG or EEG data. In a linear approximation the size of the clusters is proportional to the inverse SNR.